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Two one-dimensional spin-1 antiferromagnetic Ising models with a single-ion anisotropy 
under external magnetic field at low temperatures are exactly investigated by the transfer- 
matrix technique. The magnetization per spin (m) is obtained for the two types of models 
(denoted by model 1 and 2) as an explicit function of the magnetic field (H) and of the 
anisotropy parameter {D). Model 1 is an extension of the recently one treated by Ohanyan 
and Ananikian [Phys. Lett. A 307 (2003) 76]: we have generalized their model to the spin-1 
case and a single-ion anisotropy term have been included. In the limit of positive (or null) 
anisotropy (D > 0) and strong antiferromagnetic coupling (a = Ja/Jf > 3) the m x H 
curves are qualitatively the same as for the spin S = 1/2 case, with the presence of only one 
plateau at m/m sat — 1/3. On the other hand, for negative anisotropy (D < 0) we observe 
more plateaux (m = 1/6 and 2/3), which depend on the values of D and a. The second 
model (model 2) is the same as the one recently studied by Chen et al. [J. Mag. Mag. Mat. 
262 (2003) 258)] using Monte Carlo simulation; here, the model is treated within an exact 
transfer-matrix framework. 
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I. INTRODUCTION 

For some one-dimensional (ID) antiferromagnets at low temperatures, it has been observed a 
spin gap, which is induced by a finite magnetic field, and a plateaux structure appears in the mag- 
netization process. Experimentally, the magnetization plateaux were observed in high-field mea- 
surements of several magnetic materials such as the quasi one-dimensional compounds SrC^Os 1 , 
Y 2 BaNi0 5 2 , Ni(C 2 H 8 N 2 )2N0 2 C104 (abbreviated NENPjii, and Cu(N0 3 )22.5H 2 C£, the triangular 
antiferromagnets CeEu^, CsCuCi^ and RbFe(Mo04)2§, and the quasi two-dimensional compound, 
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with a Shastry-Sutherland lattice structure, SrCu2 (603)2-- The mechanism for the appearance of 
these magnetization plateaux in quasi one-dimensional spin chains are dimerization, frustration, 
single-ion anisotropy, periodic field and so on. 

From a general view point, Oshikawa et al.— concluded that the necessary condition for the 
magnetization plateaux in spin-S" chains is Q(S — m) =integer, where Q is the spatial periodicity 
of the magnetic ground state and m is the magnetization per site. For some range of the magnetic 
field H (i.e., H\ < H < H2), the system ceases responding to its increase and a plateau is formed in 
the magnetization versus the magnetic field curve. The values of m at which the plateaux appears 
are sensitive to small changes in the parameters of the model and are not only restricted to integer 
spin (Haldane conjecture) 1 ^. 

In the 5 = 1/2 antiferromagnetic Heisenberg model on a triangular lattice, a magnetization 
plateau was found at m/m sa t = 1/3^^^. In an S = 1/2 trimerized Heisenberg model^, the 
plateau appears at m/m sat = 1/6. Recently, plateaux at m/m sa t = 1/8 and 1/4 have been observed 
in the SrCu2(B03)2^, which has a Shastry-Sutherland lattice structure. However, irrational values 
have not been found, at least so far. Theoretically, various other models with spin S = 1/2 have 
been proposed to describe the magnetization plateaux. One of the first models was introduced by 
Hida 14 , where a Heisenberg chain was considered, with antiferromagnetically coupled ferromagnetic 
trimers (p = 3). The three-dimerized Hamiltonian proposed by Hida to describe the 3CuCl2-2 
dioxane compound is given by 

-j^trim _|_ yjint _|_ j_^Zeeman ^-q 

with 

H trim = _j F J2(S i -r i + r i -a i ), (2) 

i 

H int = J A Oi ■ S i+ i, (3) 

i 

and 

H z ~ = -fi B H {Sf + a\ + Tf) , (4) 

i 

where J a and Jp are the antiferromagnetic and ferromagnetic interactions, respectively, Sj, and 
&i are the S = 1/2 spin operators at site i, [is is the Bohr magneton and H is the magnetic field. 
Using exact diagonalization of finite systems, Hida obtained, for Jp comparable to or smaller than 
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J A, a plateau at m/m sa t = 1/3. The plateau mechanism was considered to be a purely quantum 
phenomenon, where the concepts of magnetic quasiparticles and strong quantum fluctuations are 
regarded to be of major importance for understanding the process. On the other hand, Ohanyan 
and Ananikian^ have recently studied the Hida model by using the transfer-matrix technique, 
replacing the spin operators (Sj,Tj and <7j) by Ising variables (Sf, of, rf ). It was shown that, 
for this classical model and for T = (ground state) and J a > 3Jp (strong antiferromagnetic 
coupling), a magnetization curve with plateau at m/m sa t = 1/3 is observed, indicating that the 
appearance of plateaux is not a quantum manifestation, but may be caused by the stability of 
spatially modulated spin structures. 

Another model which presents magnetization plateaux is the one-dimensional spin-1 antifer- 
romagnetic Heisenberg with single- ion anisotropy 16 . This model is described by the following 
Hamiltonian 

H = JJ2 Si ■ S m - ii B H £ Sf + L> £ {Sf? , (5) 

i i i 

where D is the single- ion anisotropy. For D = 0, the ground state is a singlet and the lowest 
excitation is a triplet (Haldane conjecture^); increasing D, the triplet splits into a higher-energy 
singlet and a lower-lying doublet, with the Haldane gap for D = 0, A(0), spliting into two gaps, as 
observed in neutron scattering of NENP 17 . The Haldane gap for general D, A(D), presents two 
different behaviors: for D > D c = J, it increases with D, while for D < D c A(D) decreases as D 
increases. 

Recently, spin S > 1 Ising antiferromagnetic chains with single-ion anisotropy have been studied 
by using classical Monte Carlo simulation 18 and it was observed the presence of 2S + 1 plateaux 
for D > 0. Essentially, these classical models are obtained replacing the spin operators (Sj) by 
Ising variables (5f) in Hamiltonian (5). From a theoretical point view, the model studied by Chen, 
et. al.— represent the ID antiferromagnetic Blume-Capel modet^, and it was observed also two 
different critical behaviors, which depend on the anisotropy parameter D (D < D c and D > D c , 
where D c = J). 

The purpose of this work is to obtain exact results for two classical models with spin S = 1 
and in the presence of a single-ion anisotropy. In Section 2 the ID models are presented and 
exactly solved by the transfer-matrix technique. The magnetization plateaux and ground-state 
phase diagrams are discussed in Section 3. Finally, the last section is devoted to conclusions. 
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II. MODELS AND FORMALISM 

The transfer-matrix technique was proposed years ago by Kramers and Wannie r 20 ' 21 , and it 
formed the basis for Onsager's solution 22 of the two-dimensional Ising model. In this section, we 
use this technique to obtain exact results for two one-dimensional models, in order to analyze the 
magnetization plateau mechanism. 

A. Model 1: Three-dimerized chain 

The first model we study is described by the following Hamiltonian: 

Hl = -j F J2(St-ri + af-Ti-aa!-S! +1 )-f, B H^2(S! + ri + a!)- 

i i 

Dj2[(S!f + (ri) 2 + (af) 2 ], (6) 

i 

where a = J a/ Jf and the spin variables Sf,T? and of can assume the values —1, 0, 1. The above 
Hamiltonian represents a nonuniform spin system in which ferromagnetic trimers composed of 
S = 1 spins (Sf,T? and of) are coupled antiferromagnetically in one dimension, in the presence 
of a magnetic field (H) and single-ion anisotropy (D). In the limit a — > (strong intratrimer 
ferromagnetic interaction), the variables S£,t* and of form a single spin £j with magnitude 3. 
Thus, the system can be approximated by a spin 5 = 3 antiferromagnetic Blume-Capel chain. 

The transfer- matrix technique is based in the calculations of the eigenvalues {Aj}, determined 
from the solution of the secular equation 

det(W x - XI) = 0, (7) 

where I is the identity matrix 3x3 and W\ the Wannier matrix, with the elements defined by 

W 1 (S,S')= exp[a(r)5 + dS 2 + b{o)S' + c{t,o)}, (8) 

cr,T=0,±l 

with 

o(r) = /3J f t + Pij-bH, (9) 
b(o) = -a/3J F o, (10) 



c(r, o) = /3J f ot + PvbH(t + o) + /3D(t 2 + o 2 ), 



(11) 
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and 

d = (3D, (12) 

where S, S' = 0,±1. 

Using properties of the matrix trace, the partition function Z = Tr(W N ) can be written as a 
sum of the iVth power of the eigenvalues {Aj} obtained from Eq.(7), i.e., 

Z = j^\f. (13) 

i=i 

In the thermodynamic limit (N — > oo), the free energy, magnetization, magnetic susceptibility 
and specific heat (per atom) are expressed in terms of maximum eigenvalue A max , respectively, as 

f = zIl nXmax , (14) 

m — — ^^ max (15) 
3A max dH 

dm _T d 1 d\ max 
X OH 3&H" V A max dH h K ' 

and 

_ 2T (9A max 2 d , 1 ^Amax , s 

C "3A max 9T + 9T l A max <9T h 1 j 

where the factor 1/3 was introduced because there are three spins in each site of the chain, and 
the maximum eigenvalue A max is given by 

Amax = -| + 27QcOs(^), (18) 

with 

A = W^l, 1) + W^O, 0) + Wi(-1, -1) = Tr(Wi), (19) 

(20) 

5 = W^(l, 0)^i (0,l) + ^i(l,-l)^i(-l,l) + Wi(-l, 0)^i (0,-1) - 

W^l, l)Wi(0, 0) - Wi(l, l)Wi(-l, -1) - Wi(0, 0)Wi(-l, -1), (21) 

<L4ff - 27C - 2A 3 
R= 54 > (22) 
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C = -det(Wi), (23) 

and 

6 = arccos(^). (24) 

Replacing the value of A max given by Eq. (18) in Eqs.(14)-(17) we obtain all the thermodynamic 
properties of a classical three-dimerized chain with single-ion anisotropy. 



B. Model 2: Antiferromagnetic Blume-Capel chain 

The second model we study is described by the following Hamiltonian: 

n 2 = jJ2s*s? +1 - n B H £ St + D^iSf?- (25) 

i i i 

which is treated also through the transfer-matrix techniqu e) 20 ' 21 . 
The Wannier matrix elements in this case are given by 

W 2 (S, S') = exp(-(3JSS' - pDS 2 + PubHS), (26) 

where the maximum eigenvalue A max is analogous to Eq.(18) with Wi(S, S') replaced by Wz^S, S') 
in Eqs.(19), (21) and (23). 

The model given by Eq. (25) has been recently treated using Monte Carlo simulation^. Using 
the transfer-matrix technique, which allows for an exact solution of the model, we can compare 
the results obtained by both procedures cited above. 



III. RESULTS AND DISCUSSION 



At low temperatures, the behavior of the magnetization for model 1 as a function of the magnetic 
field depends on the ratio of coupling constants, a = Ja/Jf, an d on the reduced anisotropy 
parameter, 5 = D/Jp. 

For S > 0, the qualitative results are the same as those obtained by obtained Ohanyan and 
Ananikian^, namely: for a > a c (5) (strong antiferromagnetic coupling) and T = a magnetization 
plateau appears at m = 1/3, for magnetic fileds in the interval H G [H Cl , H C2 ], where h Cl = 
H ci /Jp = 2.0 and h C2 = H C2 /Jp = a — 1.0. We have obtained, in this strong antiferromagnetic 
regime and positive anisotropy, that the value of the critical ratio a, a c , does not depend on 5, i.e., 
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a c (8) = 3.0. For infinity anisotropy (5 — > oo) our results reduce to the case of the three-dimerized 
Ising chain with spin 1/2. So, in this case, we obtain the same magnetization plateaux as in Ref.14. 

The ground state (T = 0), in the absence of a magnetic field (H = 0), is the antiferromagnetic 
spatially modulated structure, in which trimers of spins pointing up (Si = 1) alternate with trimers 

of spins pointing down (Si = —1) (i.e., TTTITITTT and so on), for all values of 5 > and 

a. We denote this modulated phase by (3). In the low field region (H < H C1 ), no magnetization 
is observed (m = 0). When the magnitude of the external magnetic field increases, at the critical 
value H Cl the system passes from its ground state (3) to the novel spatially modulated structure 
(3111), in which the periodic sequence of spins consists of one trimer pointing along the field, in the 
spin state Si = 1 (Ttt), and another trimer with alternating orientation of spins (|TI)- In this spin 
state (3111), m = 1/3 and increases discontinuously to the saturation value m = 1 at the second 
critical field h C2 = a — 1.0. In the above one-dimensional model, the indispensable condition for 
the appearance of the plateau at m = 1/3 is the strong antiferromagnetic coupling, characterized 
by a > 3.0. Therefore, no plateau is found for a < 3.0. 

On the other hand, when the anisotropy is negative the spin state Si = (represented by 
O) is energetically favorable, when compared to the spin states Si = 1 and —1. For certain 
values of a and 5 < 0, we can observe various types of phase transition: first, a flip of the 
central spin to the state Si = with a modulated structure (3101) (corresponding to the spin 

configuration TTT-i- O ITTTI O I ) and magnetization m = 1/6; second, a flip of the 

central spin to the state Si = 1 with a modulated structure (3111) (corresponding to the spin 

configuration TTTITITTTITI ) an d magnetization m = 1/3; third, a flip of the surface spins 

to the state Si = with a modulated structure (3010) (corresponding to the spin configuration 

ttt O T O TTT O T ••••) and magnetization m = 2/3, and, finally, the saturated state (3 2 ) 

(corresponding to the spin configuration TTTTTTTTT )> with a magnetization m = 1. 

In order to obtain the values of the critical fields h ci (transition between the modulated structure 
(3) and (3101)), h c , 2 (transition between the modulated structure (3101) and (3111)), h C3 (transition 
between the modulated structure (3111) and (3010)) and h Cs (transition between the modulated 
structure (3010) and <3 2 }), we compare the energies for the respectiv periodic sequence, finding 
the following critical fields: 
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h Cl = 2 + S 
h CQ = 2 - S 

(27) 

h C3 = a + 5 — 1 
h Cg = a — 5 — 1 

By solving Eq.(26), we obtain the critical frontiers which separate the various modulated phases, 
corresponding to magnetization plateaux at m = 1/6, 1/3 and 2/3 (see discussion in the previous 
paragraph). We find that, for < h < h c \ and h > h cs , the magnetizations are m = (disordered 
state) and m = 1 (saturated state), respectively. Note that for 5 < —2.0 there is no disordered 
state for positive field (H > 0). Depending on the values of the parameters a and S < 0, we can 
have various magnetization plateaux at m = 1/6, 1/3 and 2/3. In Figs 1, 2, and 3 the ground 
state phase diagrams for S = —1.0,-2.0, and —3.0 are depicted, respectively, where we indicate 
the various magnetization plateaux. 



IV. CONCLUSIONS 



We have treated two soluble models: model 1 (Eq. (6)) and model 2 (Eq. (25)), by using a 
transfer-matrix technique. In the thermodynamical limit (N — > oo), the partition function was ob- 
tained and ground state phase diagrams were calculated, in order to analyze the plateaux structure 
in the magnetization. In contrast to the majority of existing approaches (numerical diagonaliza- 
tion) to treat the problem of magnetization plateaux in quantum models (see, for example, Ref.13), 
the present formalism (transfer-matrix technique) is entirely based on analytical (exact) calcula- 
tions and allows for the calculation of magnetization profiles for arbitrary finite temperatures and 
values of the parameters a and 6. 
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